TENSOR ALGEBRAS, DISPLACEMENT STRUCTURE, 
AND THE SCHUR ALGORITHM 



T. CONSTANTINESCU AND J. L. JOHNSON 

Abstract. In this paper we explore the connection between tensor algebras and 
displacement structure. We describe a scattering experiment in this framework, we 
obtain a realization of the elements of the tensor algebra as tranfer maps of a certain 
class of nonstationary linear systems, and we describe a Schur type algortihm for the 
Schur elements of the tensor algebra. 



1. Introduction 

It was recently showed in || that the tensor algebras have displacement structure. 
The goal of this paper is to further develop this connection. The displacement structure 
theory, as it was initiated in the paper ||, consists in recursive factorization of matrices 
with additional structure encoded by a so-called displacement equation. Many applica- 
tions of this theory were found in applied fields, like in the study of wave propagation 
in layered media, or filtering and modeling of nonstationary processes (see || for a 
survey). On the mathematical side, there are applications to bounded interpolation in 
upper triangular algebras (including the classical setting of bounded analytic functions 
on the unit disk) and to triangular factorizations of operators. 

On the other hand, the tensor algebra is a classical concept and in H it is showed 
that there is a connection expressed by the fact that any element T of the unit ball of 
the algebra (with respect to an appropriate operator norm) satisfies the equation 

JV 

A~Y,FkAF* = GJG*, 

k=l 

where A — I — T*T and the so-called generators F^, G and J are described in the 
next section. This equation is of displacement type so we can develop a theory for the 
tensor algebra that parallels the one dimensional case. It is the main goal of this paper 
to describe several such results. 

The paper is organized as follows. In Section 2 we review the basic elements relevant 
to the tensor algebra and to the so-called noncommutative analytic Toeplitz algebra 
which is a natural extension of the tensor algebra introduced in |TT[| , and studied in 



other papers such as @, 0, JTT|, [Q. We describe an isometric isomorphic represen- 



tation of this later algebra by a certain algebra Uq-iji-) of upper triangular bounded 
operators and this represention allows us to rely on the results and ideas in the dis- 
placement structure theory. In this section we also show that a certain class of positive 

l 



definite kernels on the free semigroup with N generators, introduced in p2| , also has 
displacement structure of the type mentioned above and in Corollary 2.5 we explain 
the connection between this class of positive definite kernels and the algebra Ut(H) 
which is well-known for N = 1. In Section 3 we obtain results corresponding to so- 
called scattering experiments. Theorem 3.1 is the main result in this direction and as 
a by-product we can obtain Schur type parametrizations for Ut(H.) and the realization 
of the elements of the algebra as transfer maps of a certain class of nonstationary linear 
systems. Section 4 develops a Schur type algorithm for Ut(H). We plan to explore ap- 
proximation properties (Szego type theory) related to this algorithm in a future paper. 
This paper is part of fTj. 

2. Preliminaries 

We briefly review the construction of the tensor algebras and we explain their dis- 
placement structure. 

2.1. Tensor Algebras. We introduce some notation relevant to tensor algebras and 
we describe their representation as algebras of upper triangular operators (see, for 



instance, |10[). The associative tensor algebra T(H) generated by the complex vector 
space 7i is defined by the algebraic direct sum 

T(H) = ® k > H® k , 

where 

H m = H ® . . . ® H 



s — 



k factors 

is the /c-fold algebraic tensor product of 7i with itself. The elements of T{7i) are 
terminating sequences x = (x ,Xi, . . . ,x k ,0, . . .), where x p e Tt® p is called the pth 
homogeneous component of x (note that H®° = C, Ti® 1 = TC). The addition and 
multiplication of elements in T(H) are defined componentwise: 

{X + y) n = Xn + Vn 

and 

(xy) n = ^ x k ®y h {xQ®y n = y n <&XQ = x Q y n ). 

k+l=n 

In order to deal with displacement structure matters it is convenient to use the fact 
that T{7i) can be represented by upper triangular matrices of a special type. For a 
certain simplicity, from now on we take H = C N , N > 1. Then we denote by F(H) 
the full Fock space associated to 7i, that is, the Hilbert space 

F{H) = ® k > H m 

obtained by taking the Hilbert direct sum of the spaces 7i® k , k > 0, on which we 
consider the tensor Hilbert space structure induced by the Euclidean norm on 7i = 
C N . If {ei, e2, • • • , ejv} is the standard basis of C^, then {e^ (g> . . . e* fc | ii, . . . ,i k € 
{1,2,... ,N}} is an orthonormal basis for H® k . It is more convenient to write e CT 
instead of ® . . . ei k , where a — i\ . . . i k is viewed as a word in the unital free semigroup 
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with N generators 1, . . . ,N. The empty word is the identity element of F^. The 
length of the word a is denoted by \a\. We also use to write the elements of F^ 
in lexicographic order. The space 7i® fc+1 can be identified with the direct sum of iV 
copies of TL® k . More precisely, we notice that {{Sji^a)f=i I j £ {1 } • • • ;^} ; W\ — k} is 
an orthonormal basis for H® k © . . . © H® k j = (U m ) m (8ji is the Kronecker symbol). 

N terms 

We deduce that the mapping 

^k{{{5jie a )f =x )) = e ja , j = l,... , N, 
extends to a unitary operator from (}{® k )® N onto 'H® k+1 and the formula 

gives a unitary operator from Tik = 7i® Nk 1 onto 7i® k (the notation T® 1 means the 
direct sum of I copies of the operator T). We finally obtain a unitary operator ip : 
©£L 7"4 - * T(J-C) by the formula if) = ffi^Lo^fc, where ip is the identity on C = 7i - 
Each Tiki k > 2, is also identified with the direct sum of N copies of Tlk-i, 

H k = H® Nh - 2 @...@H® Nk -* 
^ j 

V 

iV terms 

but the explicit mention of the identification can be omitted in this case. 

An operator T G £(©£L 7ifc) has a matrix representation T = [Tij}°° =0 with T^- G 
C(Hj,Hi). We say that T is upper triangular provided that Ty = for z > j. We 
denote by U^-(TC) the set of upper triangular operators T = [Tij]°° =Q in £(©£i 7ifc) 
with the property that for i < j and i, j > 1, 

(2.1) r„ = _! 

and To-,- = for j > k and some fe. We notice that the entries To-,-, j > 0, determine 
the matrix T. Also, a simple calculation shows that U^-(TC) is an associative algebra. 

We define an algebra isomorphism 

(2.2) $ : T{H) -> Z-#(W) 

as follows: for x G T(H), x = (x ,Xi, . . .), there are complex numbers c CT , cr G F^, 
such that 

(2.3) Xj = ^ ^ Co-Co-. 

For j > 0, T j is the row matrix [co-] [0-1=7. Then T j = for sufficiently large j and we 
can define T G C(®^L Ttk) by using Q2.1|). Set $(x) = T. Clearly $ is a bijection. 
In order to see that $ is also an algebra morphism it is convenient to introduce the 
following operators. First, the operator C~, j = 1, . . . , N, is an element of Uq-{7i) 

defined by the formula: Cj = $(0, ej, 0, . . . ) = [T^.]^._ , where Tj fc = for k ^ 1 and 

Tq X = Cj. Then Cj~ = (C~)*, j = 1, . . . , N. If we use the notation C~ to denote the 
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operator C ix . . . C i , where a — %\ . . . ik, then we deduce that each T G U^-{H) has a 
representation 

(2.4) T=£c CT C- 

Now the required properties of $ follow from the easily verifiable fact that 

$((0,e„0,...)®(0,e fc ,0,...)) = CrC,7. 

It appears to be quite natural to consider the algebra Ur(H) of all those opera- 
tors in £(©£l 7ifc) that satisfy the condition ( [2.1| ). It turns out that Ur{Ti) is iso- 
metrically isomorphic to the so-called non- commutative analytic Toeplitz algebra, re- 
cently studied in papers like |TTJ, 0, ||. We also denote by S{7i) the Schur class 
of all contractions in Ur^H). We can extend this setting as folows. Let £\ and £2 
be Hilbert spaces, then V(q-(7~C, Si, £2) denotes the set of the operators T — P~y?°_Q 
in C(®^ =0 (H k £1), ®%L {H k ® £ 2 )) obeying flp). Instead of U T {H,E,E) we write 
Ur(H,£)- Also S(H,£i,£2) denotes the set of the contractions in Ut( / H,£i,£2)- We 
notice that U r {H) = U r {H, C, C). 

2.2. Displacement Structure. The systematic study of the displacement structure 
was initiated in || , and since then manifold applications have been found in the fields 
of mathematics and electrical engineering (for a recent surveys, see 0). Two main 
themes of this theory are about scattering experiments and recursive factorizations. A 
matrix R has displacement structure with respect to the generators F and G provided 
that 

(2.5) A - FAF* = GJ pq G*, 

where J pq = I p @ —I q is a signature matrix that specifies the displacement inertia of A. 
In many applications, r = p + q is much smaller than the size of A. From now on we 
will write J instead of J\\. There is also a version of this theory allowing A to depend 
on a (usually integer) paramater, see || or ||. A major result concerning matrices 
with displacement structure is that, under suitable assumptions on the generators, the 
succesive Schur complements of A inherit a similar structure. This allows to write the 
Gaussian elimination for A only in terms of generators (see ||, 0). 

We now describe a connection between displacement structure and Ur(H)- Let 
Cfc = [Z^]°° =0 , k = 1, . . . ,N, be isometries defined by the formulae: = for 
i 7^ j + 1 and is a block-column matrix consisting of blocks of dimension 

dim7ij, all of them zero except for the kth block which is 1%. These operators are 
closely related to the operators C£, in the sense that there is a unitary operator u such 
that Cfc = uC^u*, k = 1, . . . , N. The following result was noticed in [j5|. 

Theorem 2.1. Let T = [Ty]^ G S(H) and A = I - T*T. Then 

N 

(2.6) A - CkAC* k = GJG*, 

k=l 
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where the generator G is given by the formula: 



G 



1 T ( 
T r 



oo 



01 



The equation (|2.6|) can be rewritten in the form of a time variant displacement 
equation by setting F = [C\, . . . Cn}- We now discuss some additional examples. 

2.2.1. Positive definite kernels on F^-. A mapping K : x F^ — > C is called a positive 
definite kernel on F^ provided that 



2j K{a h o-j)\j\i > 

for each positive integer k and each choice of words o~\ , . . . , in F^ and complex 
numbers Ai, . . . , A*,. Without loss of generality, we assume K(0, 0) = 1. The semigroup 
F^ acts on itself by juxtaposition and we assume that the positive definite kernel K is 
invariant under this action, that is 

K(ra,ra') = K(a,a'), r,a,a' <E F+. 

By the invariant Kolmogorov decomposition theorem, see e.g., [10|, Ch. II, there exists 



an isometric representation u of F^ on a Hilbert space /C and a mapping v : F^ — > /C 
such that 

K(ct,t) = (v(t),v(o-)) k , 
u(t)v(o~) = v(to) 

for all a,r E F^, and the set {v(a) \ a G F^} is total in /C. The mapping v is unique 
up to unitary equivalence and u is uniquely determined by v. Also, the previous 
representation of K implies, as noticed in |T2|] , that u(l), . . . , u(N) are isometries with 
orthogonal ranges if and only if K(o~, r) = for all pairs (a, r) with the property that 

ar or r 



there is no a E F^ such that a 



aa. In this case, K is called in 12 



multi-Toeplitz kernel and we will use the same terminology. We can define 



(2.7) 



By induction on the length of words and using properties of the lexicographic order 
on words we check that if K is a positive definite multi-Toeplitz kernel, then for all 
i, 3 > 1, 



(2- 



Si 



This relation is an analogue of (|2.1|) , therefore we expect that positive definite multi- 
Toeplitz kernels have displacement structure. 



Theorem 2.2. Let S 



oo 



6e t/ie array of operators associated to a positive 



definite multi-Toeplitz kernel on W%. Th 



en 



N 



(2.9) 



c kSC* k = GJG* 



k=l 



where the generator G is given by the formula: 

1 

G = $01 "oi 



The equality above can be interpreted in the sense of an equality of unbounded oper- 
ators with dense domain made of the terminating sequences in (B^qHi- 

Proof. We write S = L + Q, where L = [Ly]°° =0 and Q = [Qij}°° =0 are the lower and, 
respectively, the upper parts of S, such that = Sy for i > j and = for i < j, 
while Qij = Sij for i < j and = for i > j. Since, formally (that is, computing the 
entries of the arrays that are involved according to the usual rule of multiplication of 
operators), LCk = C k L and QCl = CIQ for all k = 1, . . . , N, we deduce that 

S-TZ±C k SCi = L + Q-^ =1 C k (L + Q)C* k 



L + Q - Y^k=i CkLCl - Ylk=i CkQC t 



N 



L + Q - Y,k=i LC k C* k - Ylk=i CkClQ. 



N 



But / = J2k=i CkCl = Pn , the orthogonal projection on Hq, so we can deduce that 



~~ l^k=l U fc OL/ fc 



L + Q-L(I-P Ho )-(I-P Ho )Q 
LP-Hq + Pu Q — GJG* . 



□ 



We now consider truncations of a positive definite multi-Toeplitz kernel on F^. If K 
is such a kernel and are associated to K by the formula ( |2.7| ), then we define K n = 
[*%]ij=o- These are positive operators on ©[L 7iz and their entries satisfy ( |2.8| ). We 
also introduce the operators C k)U = P®n_ o ^ z C) 
of a displacement equation. 



=0 Ti.i- We obtain another example 



Corollary 2.3. The positive operator K n satisfies the displacement equation 



N 



(2.10) 



K n ~ C kin K n C kn — GJG* , 



k=l 



where the generator G is given by the formula: 



G 



q* C* 

°01 °01 



q* 
J 0n 



q* 
J 0n 



Proof. This is an immediate consequence of Theorem E72 



□ 



It turns out that the operators satisfying ( |2.10| ) have a certain explicit description. 
For x G ffi" =0 7iz we denote by U (x) the upper triangular operator on ©" =0 7i; satisfying 
and whose first row is x. 



Theorem 2.4. A positive operator A on ®f =Q 7ii satisfies ( 2.1 0|) if and only if A = 
U(x)*U(x) — U(y)*U(y) for two vectors x,y G ©™ =0 7-^ with the property that there is 
another vector z G (B^qHi such that U(y) = U(z)U(x) and U{z) is a contraction. 



Proof. If A = U(x)*U(x) - U(y)*U(y) then we use the fact that for any x G ®[Lo^ 



and k = 1, . . . , N, C n ^U(x)* = U(x)*C n ^, and deduce ( |2 . 1 U| ) as in the proof of Theo- 
rem |2]2|. Conversely, if A satisfies ( |2.10| ), then we use the fact that C CTjn = for \o~\ > n 
in order to deduce that 

n 

A — GJG -\- ^ C an GJG C a n . 
M=i 

Hence, if G = [x*,y*] for some x,y G ©™ =0 7^, then A = U(x)*U(x) - U(y)*U(y). Since 
A is positive, it follows that there exists a contraction T such that U(y) = TU(x). 
We notice that we can assume, without loss of generality, that Xq ^ and then we 
deduce that C nt kT* = T*C Ut k for k = 1, . . . ,N, which implies that T — U(z) for some 



z G 



□ 



We mention a simple but useful consequence of this result. Let K be a positive 



definite multi-Toeplitz kernel on with K($, 0) 
operators associated to K by (|2.7|). Let 



1 and let SV,-, i,j > 0, be the 





' 1 " 




' " 




q* 
°0l 




q* 
°01 


•En 




, Vn = 






q* 

I J 0n J 




c* 

L J 0n J 



then K n = U(x n )*U(x n ) — U(y n )*U(y n ) and by the previous result there exists z n G 
©™ =0 7"^ such that U(y n ) = U(z n )U(x n ). In fact, z n is uniquely determined and also 



P+ 



n+l, 



.qTCi = U(z n ) for all n > 1. It follows that there exists a unique 



T{K) G S(H) such that P &LqHi T \ ®? =0 Ht ■ 
the set of those elements T G S{7i) with T( 



= 0. 



and T(K) 00 = 0. Denote by S (H) 



Corollary 2.5. The mapping K — » T(K) gives a one-to-one correspondence between 
the set of positive definite multi-Toeplitz kernels on iwit/i X(0,0) = 1 and So(H.). 
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2.2.2. Pick kernels. This type of kernels were recently introduced in ||, |]I2"|| , and they 
were shown to have a certain universality property with respect to a Nevanlinna-Pick 
problem, M. 

Let B N be the open unit ball in C^. For distinct points Ai, . . . , Xl in B N and 
complex numbers b\, . . . , bL, we define the Pick kernel by the formula 



R 



1 I'jhi ' lL 



1-<A,,A Z ) 



j,/=i 



where (A 3 -, A;) denotes the Euclidean inner product in C^. 
We can define 

Fj = diag(Ap ) )f=i, j = l,...,N, 

the diagonal matrix with the diagonal made of the jth components of the points Ai, 
. . . , Xl. Also, we set 

i h 

1 b 2 



G 



1 b 



L . 



We can check by direct computation that 



N 



(2-11) 



R — ^ FkRFfc = GJG*. 



k=l 



It was shown in how to use this equation in order to solve some multidimensional 
Nevanlinna-Pick type problems as those in [Ef, fT2||. 



3. Scattering Experiments 

One of the main results in displacement structure theory refers to the possibility of 
associating a scattering operator to the data given by the generators of the equation 
( p.5|) . In this section we prove a similar result for displacement equations of the type of 
( p. HQ and we explain the connection with a special class of time variant linear systems. 



Let Si, S 2 , and Q be Hilbert spaces. We consider a displacement equation of the 
form 



N 



(3.1) A - V F k AF* = GJG* 



k=l 



where F k G C(Q), k = 1, . . . ,N, are given contractions on the Hilbert space Q. Also 
G = [ U V 1 G C(Si © S2,Q) and J = Ie 1 © Ie 2 . The wave operators associated 
to ( |3.1| ) are introduced by the formulae: = [Uk]^L , V<£> = [Vk]^L , where Uk = 
[F a U] H=k :H k ®Si^g,V k = [F a V} W \ =k :H k ®S 2 ^Q. We will assume that both 



Uoo and Voo are bounded and also that lim^oo J2\ a \=k \\F£g\\ = for all g G Q. Under 
these assumptions we deduce that (|3.1|) has a unique solution given by 

(3.2) 



A = U^Uo, - CKc- 



The next result was already noticed in ||, but here we present a different approach 
based on system theoretic ideas. 

Theorem 3.1. The solution ( |3.2|) of the displacement equation ( |3.1|) is positive if and 
only if there exists T G <S(7i, £1,82) such that = TU^. 

Proof. Assume A = U^Uoo - V^Km > and let A = LL* be a factorization of A with 
L G (?) for some Hilbert space T . From (|3TT| ) we deduce that 

N 



LL* + VV* = F k LL*F* + UU* . 



k=i 



In matrix form, 



(3.3) 



[L V] 



L* 
V* 



[ F\L ... F N L U ] 



L*F* 



L*F i 
U* 



N 



Define A* = [ L V ] and B* = [ F 1 L ... F N L U ], then we deduce from 

that there exists an unitary operator 9 G £(1Z(B),1Z(A)) such that A = 9 B. It follows 
that there exist Hilbert spaces TZi, 1Z 2 , and an unitary extension 9 G £(jF ffiAr © £4 © 
TZi,^ © £2 © ^2) of 6>o, hence this extension satisfies the relation 



(3.4) 



.4 
0% 
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B 

0% 



Let i G {1,2,3}, j G {1,2,... ,N + 2}, be the matrix coefficients of It is 
convenient to rename some of these coefficients. Thus, we set 



9 lk , k 



Y k = 9 2 k, k 
From (I3T1) we deduce that 



l,...iV, 



N, 



Z — 6i,n+i, 

W = &2,N+1- 



N 



L* = X k L*F^ + ZU* 



and 



k=i 



N 



V* = J2Y k L*F* + WV*. 



k=l 



By induction we deduce that 

N n 

(3.5) v* = wu* + Y kXaZu*F* a + QtL*f;, 

k=\ |cr|=0 \r\=n+l 

where Q T are monomials of length \t\ in the variables X\, ... ,X^, Yi, . . . , Yjy Since 
9 is unitary it follows that all Q T are contractions. 

We define T 00 = W and for j > 0, 

Then we define Ty, i > 0, j > i, by the formula and = for i > j. We show 
that T = [Tij]fj =0 belongs to S(TC, 8±, £ 2 )- To that end, we introduce the notation: for 
n > 0, 



X 1 (n) =[X 1 ... X N ] eJV , ^(n) = [ Fx ... ^ 

X 2 {n) = Z® N \ Y 2 (n) = W® N " . 

We check by induction that T 00 = Y 2 {0), T 01 = Y 1 {0)X 2 (1), and for j > 2 

r oi = y 1 (o)x 1 (i)...x 1 o--i)x 2 (j). 

Consequently, T* is the transfer map of the linear time variant system 
x(n + l) = Xi(n)*x(n) +Yi(n)*u(n) 



y(n) = X 2 {n)*x{n) +Y 2 (n)*u(n), n>0. 

, n > 0, is a contraction, it follows that T* is a con- 



Xi(n) X 2 (n) 
y 2 (n) 

traction, hence T is a contraction and T G E\,E-i). Also, since lim^oo X]|<r|=fc 11-^ d\ 

for all g & Q, we deduce from (|3.5|) that = TUoq. □ 



Since each matrix 



This result is a generalization of Theorem It can be used to solve interpolation 
and moment problems in several non commuting variables as suggested in |J. The 
proof given here emphasizes the fact that this type of problems appears as an interesting 
particular case of similar questions in time variant setting (as described, for instance, 
in H). Next we develop this idea by giving a parametrization of the elements in 
S(H, £i,£ 2 ). 

3.1. Schur Parameters. We introduce a Schur type parametrization for S(7i, Si,S 2 ) 
and then give a number of applications. A first instance of such a result might be con- 



sidered Euler's description of SO(3). Schur's classical result in |13[ gives a parametriza- 
tion of the contractive holomorphic functions on the unit disk and Szego's theory of 
orthogonal polynomials provides an alternative to this result for probability measures 
on the unit circle. 



Here we provide a refinement of a parametrization of S(H, £\, £ 2 ) given in [pjj . We 
can deal with a slightly more general situation than S{7i). For a contraction C we 
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denote by T>c its defect space, that is the closure of the range of the defect operator 
Dc = (I — C^C) 1 / 2 . Let 11(81,82) be the family of sequences {7 CT } CTgF + with the 
properties: 70 G £(8i,8 2 ) and for \a\ > 0, 

(3.6) lka eC(v la ,v lL J, 

where k = 1,...N and ka — 1 denotes the predecessor of ka with respect to the 
lexicographic order. 

Theorem 3.2. There exists a one-to-one correspondence between S(7i,8i,82) and 

n(8 h 8 2 ). 

Proof. A simple proof can be obtained as an application of the Schur parametrization 
of upper triangular contractions in Let T G S(H,8i,8 2 ) and for each j > 0, T 0j 
is the row matrix [co-]^^.,-. We use Theorem 2.2.1 in M in order to associate to T a 
family T of Schur parameters. The Schur parameters of the contraction 

. 



Too 




Tox 



C0 





Cl 
C0 



CN 









c . 



rSfand r 01 eC(v® N ,V- 



are given by r 00 = 70 = C0, r n - * 00 - ui v- -v- 70 > -7 > 

For the next step we associate Schur parameters to the contraction 



We notice that 





Too 


Tqi T02 









T\\ T12 












T22 






T 12 " 




-'oo 


1 01 





T22 







rp®N 
1 11 



By the formula (2.2.3) in Q, the entries of the matrix 



are polynomi- 



al 1 T12 
T22 

als in the Schur parameters, their adjoints and their defect operators. Consequently, 
the contraction associated to the Schur parameters r® , r®i and r® is precisely 



rp®N 

-'oo 




rp®N 

i 01 

rp®N 
1 11 



. From the uniqueness of the Schur parameters, we deduce that the 



Schur parameters of 



are precisely V 



®N -p07V 



(l0 , Tq" and r® . By induction, we 



T\\ T\2 
T22 _ 

deduce that we can associate to T the family of Schur parameters V = {Tij}i<j such 
that each satisfies ( |2.1|) . We notice that Toj ■ T>® ._ l — > T>r* _ 1 for j > 1- Also, 
each r j is a row contraction, T j = [y®]\ a \=j. By Proposition 1.4.2 in |3j], 

r i = [ 71 D ltl2 ... D lt ... d Yn _ i1n ] , 



where 7^ G £(X> 70 ,X> 



k 



N . Also, there exist unitary operators U\ : 



V 



r i 



{D lk and : "Dp* — * £7* • Next, we have that T 2 : T> 



r i 



Dr» and 

1 m 



11 



we can define the row contraction r 2 = [7°][<r|=2 from (©^ 1 D 7fe ) eAr into T> 7 ^ by the 
formula f 2 = viT 2(uf N )* . As above, 



02 



[ 7n ^712 • • • D lti . . . D^ NN _ xlNN ] 



where jkr G ), k — 1, . . . , N, \r\ — 1. We proceed by induction in order 

to construct the whole family of Schur parameters {la} ae ^+ with 70 G £(81,82) and 
obeying (|3.6|). □ 



The parameters {r 0j }j> are the Schur parameters associated to T in |TT|. The above 
proof explain their construction within the framework of [|J. The new parameters 7 CT 
are more like "true" Schur parameters in the sense that they are contractions subject 
to the only restrictions given by the compatibility relations ( |3.6| ). 

A similar parametrization can be obtained for the set of positive multi-Toeplitz 
kernels. Denote by IIo^i,^) the set of families {-j a } ae¥ + in 11(81,82) with 70 = 0. 

Theorem 3.3. There exists a one-to-one correspondence between the set of 8 -valued 
positive multi-Toeplitz kernels with K($, 0) = 1 and IIo(£). 



Proof. This is a consequence of Theorem |3.2| and Corollary |2.5| . A direct proof can be 
obtained as an application of Theorem 1.5.3 in H. The parametrization given here is 



a refinement of one given in [12 . □ 




assume that for all n > 0, the operator 



is a coisometry. The evolution 



An application gives another representation of the class S(TL, 8\, 82). Motivated by 
the proof of Theorem [TT] we consider the class BSn(8\,82, T) of systems 

A(n)x(n) + B(n)u(n) 

C(n)x(n) + D(n)u(n), n > 0, 

where 8 U 8 2 , T are Hilbert spaces, A(0) G L(T , T® N ), B(0) G C(8 2 , F eN ), C(0) G 
C(F,£i), D G C(8 2 ,8 X ), and A(n) = A(n - 1)® N , B(n) = B(n - 1)® N , C(n) = 
C(n - 1)® N , D(n) = D(n - 1)® N (these relations are all of the type of (pTT|)). Also we 

" A(n) B(n) ' 
C(n) D(n) 

of the system S can be encoded by the transfer map that shows how the infinite 
input sequence (u(0),u(l), . . .) is transformed by the system into an output sequence 
(y(0),y(l), . . .) provided that x(0) = 0. It is easy to see that the transfer map Ts is 
a lower triangular contraction with the matrix coefficients T^(n,n) = D(n), Tj](n + 
l,n) = C(n+l)B(n) and for j > 1, T^(n+j,n) = C(n+j)A(n+j-l) . . . A(n+l)B(n). 

Theorem 3.4. T belongs to S(7i, 81,82) if and only if T* is the transfer map of a 
system in BSn(8i,82,T) for some Hilbert space T. 

Proof. We use the construction from the proof of Theorem 2.3.3 in [|J. Thus, let 
{Tij}i<j be the Schur parameters associated to T in the proof of Theorem |3.2| . For 

12 



each i > we consider the row contraction 

^ = [ r« D r *.r i>i+1 D T *p Tlx+ v i)i+2 ... ] . 

Let Vi be the isometry associated to Lj by the formula 1.6.3 in [|J, that is, 

— r*jr i>i+ i — r* i z? r *. +i r i)i+2 . . . 
(3.7) v l= o D r . i+1 -r* l+1 r M+2 ... 

o o D r ^ +2 



Tu B{i)* 
C(i)* A(i)* 



V is defined on the space £ 



)%L V Vok and let 

N 



%L T>®^ with values in in £ 



)r=o^r 0fc - We set 



- Ok 



be the unitary operator defined 



by the formula </>(©|=i ®T=o u ki) 



x =1 Uki- Then we define 



C*(i) 



B*(i) 

A*(i) 



7fc = ^1 = 

r 



B*(i)(j)® Ni 
)*A*(i) 



By Theorem 2.3.3 in [§, = C(i + l)B(i) = C(i + l)fl(i) and for j > i, 

(3.8) 7;* = C(j)A(j -l)...A(i+ l)B(z) = C(j)A(j -l)...A(i+ l)B(i). 

Since the elements of the family T satisfy (|2.1|) and the matrix coefficients of Vi are 
monomials in T^-, their adjoints and defect operators, it follows that the system as- 
sociated to the families {A(n)} n > , {B(n)} n > , {C f (n)} n > and {D(n)} n > belongs to 
BS(£, T) and T* is the transfer map of this system. □ 



4. Recursive Constructions 



As noticed in the proof of Theorem |3.4| , the Schur parameters introduced in Theo- 
rem 1333 are basically given by the parametrization of the unitary group. A convenient 



way to encode the dependence of the matrix coefficients of T on the Schur paramaters 



is given by a continued fraction algorithm introduce for iV = 1 by Schur in [13]. Here 
we describe such an algorithm for N > 1. Of course, the main dificulty consists in 
keeping all the contructions within the tensor algebra. Due to the fact that there is 
a standard filtration on the tensor algebra, the contructions do not go through as in 
the one dimensional case, still we can obtain a sort of graded Schur type algorithm. 
For the algebra of upper triangular matrices, a Schur type algorithm was known (as a 
consequence of a similar algorithm for the spectral factorization of a family of matri- 
ces with displacement structure, see B or H for a detailed discussion). If we use it 



directly to our situation, it will produce the paramaters {7 CT } crg F+ of Theorem |372| , but 
the objects produced at each step are no longer in the tensor algebra. In this section 
we show how to remedy this situation. 

We first notice another useful representation of a transfer map (we will assume in 
this section that £\ — £\ — C). It is convenient to extend the symbol map Q2.2|) . Thus, 
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for T £ C((B^L Ti.k <8> £i,©j£Lo^fc ® ^2) an upper triangular operator satisfying ( pTT 
we write 



Now, if T s is the transfer map of a system £ in BSn(F) = BS^(C, C, JF), then 
T* = $(D(0)*,0,...)+$(0,B(0)*,0,...) 

xci-^o,^)*^,...))- 1 ^^)*^,...), 

provided that ||$(0, A(0)*, 0, . . . )|| < 1. We notice that this is the case if and only if 



||A(0)|| < 1. Let A(0)* = [^(0)*]f =1 , 5(0)* = [B,-(0)*]? =1 . Then 



iiV 



A? 



*(o,A(o^o,..o = £*(A i (o^o,...)c; 



and 



N 



so that 



(4.1) 



$(0, B(0)*, o,...) = J2 HBj(oy, 0, . . . )Cr, 

3=1 

T* = $( J D(o)*,o,...) + (Ef=i^(or,o ) ...)C7) 

x(I - Ef=i *(^(0)*, 0, . . . )C7)-i$(C(0)*, 0, . . . 



provided that ||A(0)|| < 1. Next, we prove a refinement of (|4.1|) which is the key step 
of a Schur type algorithm. We introduce the translation operator given by the matrix 



S 





1 

1 



Theorem 4.1. IfT belongs to S(H) then there exists an unique 2\ £ S(H,T>^ ,T) T * o ) 
such that 

T = <S>(T O0 ,0,---) + HD T j o ,0,---)TiS 

x(/-$(-T * ,0,...)T 1 5$( J D Too ,0,...). 

Proof. Let {Tij}i<j be the Schur parameters associated to T as in the proof of Theo- 
rem |3.2| . Then | Too I < 1- If I Too I = 1, then we must have = for j > i, so that 
T = $(roo, 0, . . . ) and T\ — . . . — — 0. Assume |r o| < 1. Then T* is the transfer 
map of the system associated to the operators Vi defined by (|3.7|) . Thus, we have 
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Va 



r 00 s(oy 
c(oy A(oy 



Too 


d t * r 01 

1 00 ul 


Dr* Dr* l^n? . . . 

1 00 1 01 uz 


D T 00 


_r o roi 


-^00 D T* 01 ^02 ■■■ 





Dr 01 


-^^02 








Dr 02 



Too 


Dr* 

1 00 


... " 




" / 








D r 0Q 


■p* 
1 00 


... 







r i 










/ 







^r 01 





Dr* 



— r* n 

1 00 u 

o / 



Too 


Too 

^r 00 
C"(0) 



JO 

o r i B'(oy 
o c(oy A'(oy 



D n?oi D r * o B'(oy 
v* 00 B'(oy 
A'(oy 



r oo r oi 



We assume first that ||A(0)|| < 1. Then formula (|4.1|) gives 



T = $(r 00 ,0,...) + $(0, J B(0)*,0,...) 

x(I-$(0 1 A{oy,0,...))- 1 *(C(oy,0,...) 

= $(r 00 ,o,...) + [ $(o,£> rSo roi,o,...) $(o,D r * o B'(oy,o,---) ] 



i - $(o, -r* r 01 , o, . . . ) -$(o, -r* B'(oy, o, . . . 
-$(o, c'(oy, o, . . . ) / - $(o, A'{oy, o, . . . ) 



*(D roo ,0, 




Since ||A(0)|| < 1, we also have ||A'(0)|| < 1 so that we can use a well-known formula 
for the inversion of a 2 x 2 matrix in order to deduce that 



$(r 00 , 0, . . . ) + <S>(D r * o , 0, . . . )($(0, r 00 , 0, . . . ) A- 1 

+$(o, B'(oy, o, ...)(/ - $(o, A'{oy, o, . . . j)- 1 



x$(0,C"(0)* ! 0,...)A- 1 )$( J D roo ,0,...; 
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where 

a = /-$(o,-r* r 01 ,o,...)-<f(o,-r* OJ B'(o)*,o,...) 

x(J - $(0, A'(0)*,0, . . .))- 1 $(0,C"(0)*,0, . . .) 

= /-$(o,-rs ,o,...)($(o,-r 00 ,o,...) 

+$(o, B'(oy, o, ...)(/ - $(o, A'(oy, o, . . . ))- 1 $(o, c"(o)*, o, . . . )) 

Since ||$(0, — Tq , 0, . . . ) || < 1, A is invertible and the previous calculation makes sense. 
Also, we define 



r = $(o, r i, o, . . . ) + $(o, B'(oy, o, ...)(/ - m A'(oy, o, . . . ))- 1 $(o, c'(oy, o, . . . ). 

Since 



r i B'(oy 
c(oy A'(oy 

that R e S (H). Actual 



is a contraction (in fact, it is an isometry), we can easily check 

y, we can notice that R is obtained from the Schur parameters 
{Tij}i>j, where Fa = and fy = for j > i. We obtain that 

(4.2) T = $(r 00 , 0, . . . ) + $(D rSo , 0, . . . )R(I - $(-r* , 0, . . . ) J R)- 1 $( J D roo , 0, . . . ). 

Since ||$(— Tq , 0, . . . )|| < 1, the previous formula makes sense without our assumption 
that 11^4(0)11 < 1. In fact, an approximation argument shows that (4.2) holds without 
that assumption. Thus, let p G (0, 1) and consider T* the transfer map of the system 
{pA(n),B(n),C(n),D(n)} (not necessarely in BSn(H)). All the previous calculations 
go through and we obtain a contraction R p such that ( (4.2[ ) holds for T p . Letting p — > 1 
we obtain (gg). □ 

This result leads to the following Schur type algorithm for a T e S(TC). 

T = T, 7„ = T„ j00 , n > 0, 

^ T n+1 = the unique solution of the equation 
T n = $( 7n ,0,...) 

+$(^,0, . . . )T n+1 s® Nn (i - $(- 7 ;, 0, . . . )T„ +1 ^' l )- 1 $( J D 7n , 0, . . . ). 

The algorithm generates a sequence of elements T n e S(T-C,T> t bn ,Vt* ) and it is 

n — 1 n 1 

easily seen that 7o = r 0n , n > 0. Therefore {7 n }n>o uniquely determine T and the 
algorithm encodes the dependence of T on these parameters. 
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